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The axial next-nearest-neighbor Ising model is studied in two dimensions at finite temper- 
ature using the density matrix renormalization group. The model exhibits phase transition of 
the second-order between the antiphase in low temperature and the modulated phase in high 
temperature. Observing the domain wall free energy, we confirm that the modulation period 
in high-temperature side is well explained by the free-fermion picture. 
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1. introduction 

Periodically modulated structures may occur in a 
wide range of physical systems. As examples of such 
systems, La6CagCu2404i and Ca2Y2Cu50io are well 
known, 1 ' 2 where spins of the copper atoms interact 
ferromagnetically between the neighboring sites along 
the CuC>2 chains and antifcrromagnetically between the 
next-nearest-neighboring ones. A phase transition of 
commensurate-incommensurate type was observed in 
these systems. Another example is cerium antimonide 
(CeSb) 3 which has a nontrivial phase diagram and which 
shows modulated spin patterns with various periodic- 
ities. In some ferroelectric materials, such as NaNO,3, 
the modulated phases are present between the ferro- 
electric low-temperature state and the paraelectric high- 
temperature one. 4 ' 5 

Physical properties of magnetically modulated struc- 
tures can be described by simplified models with com- 
peting interactions. One of the simplest examples is 
the so-called axial next-nearest-neighbor Ising (ANNNI) 
model, which contains ferromagnetic coupling J 1 be- 
tween nearest-neighbor spin pairs and antiferromagnetic 
one J 2 between next-nearest-neighbor spin pairs in a pre- 
ferred direction. 6 Several analytical methods have been 
developed to study the phase diagram of the ANNNI 
model in two dimensions. For instance, the free-fermion 
approximation treats domain walls running along the 
chain direction. 7 ' 8 The Muller-Hartmann-Zittartz ap- 
proach assumes existence of the domain wall in the per- 
pendicular direction to the axial one. 9 A detailed survey 
of earlier works on this topic has been reviewed by Selke. 6 
Recent progress can be found in Refs. [10-12]. 

In this paper we focus on the two-dimensional (2D) 
ANNNI model, which is described by the Hamiltonian 

H = - Jl ^2 a i,j( a i+l,j + ~ ^2 ^ 0'i,jO'i+2,j (1) 

hj hj 

on a square lattice, where the index i specifies the po- 
sition along the axial direction. The Ising spins <7ij = 
t or | interact ferromagnetically (Ji > 0) between the 
nearest neighbors and antiferromagnetically ( Ji < 0) be- 
tween the next- nearest neighbors. The ratio between the 
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Fig. 1. The ordered phases of the 2D ANNNI model. 



coupling constants k = —J2/J1 is commonly used for the 
measure of the frustration. It is widely accepted that in 
the low temperature region the model shows a ferromag- 
netic structure when k < 0.5, and when k is larger than 
0.5, the so-called antiphase structure {•■•TTI4TT'''}i s 
realized. 6_10 ' 15 ' 19 ' 20 Figure 1 shows the location of these 
ordered phases. It has been confirmed that these ordered 
phases are bordered by the second order phase transition 
lines. 

There is an argument about the presence of incom- 
mensurate (IC) phase in the highly frustrated region, 
which is specified by the condition k > 0.5. Though a 
wide area of the IC phase is expected by the mean-field 
theory, the Monte Carlo (MC) simulation by Sato and 
Matsubara suggests that the region of the IC phase is 
very small. 15 Recently, Shirahata and Nakamura per- 
formed an extensive calculation by use of the non- 
equilibrium relaxation method. 10 Assuming the pres- 
ence of the Berezinskii-Kosterlitz-Thouless (BKT) tran- 
sition 13, 14 they estimated two critical temperatures bor- 
dering the IC phase. What they found is that these two 
transition temperatures are almost identical. They spec- 
ulated that successive phase transitions may occur within 
an infinitesimally narrow temperature region. Table 1 
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summarizes these theoretical and numerical estimates 
of the phase transition temperatures at k — 0.6, where 
the T c represents the upper border of the antiphase, and 
where T" is the lower border of the paramagnetic phase. 
(The IC phase is present if T' is larger than T c .) 

The aim of our study is to obtain the precise mod- 
ulation period of the local magnetization and its decay 
factor in the parameter region where the presence of IC 
phase has been discussed. For this purpose we employ 
the density matrix renormalization group (DMRG) 16-18 
method, and carry out a scaling analysis on the domain- 
wall free energy. As shown in the following, we confirm 
that the modulation period is well explained by the free- 
fcrmion picture. 

2. Application of DMRG 

We consider the 2D ANNNI model on the square 
lattice of the size L x oo. The transfer matrix 



of this system T L [a'\cr] connects two 



adjacent spin 
rows [a'] = {a lj ,a 2j ,...,a L j} and [a] = {o^-i, 
<72,j-i> ■ • • > a L,j-i}i wncre index i runs from 1 to L to- 
ward the axial direction. For simplicity, we drop out the 
indices j and j — 1 from the Ising spin variables in the 
following, and write them as [a'] = {a' x ,a' 2 , . . . ,a' L } and 
[a] = {a 11 a 2 , . . . ,a L }. Without loss of generality, the 
transfer matrix can be written as the product of the over- 
lapped local weights 

L-2 

t l W>] = II ^K'^+i^+2k^ + i^+2) - (2) 

1=1 

where W(a' i a' i+1 a' i+2 \a i a i+1 cr i+2 ) is the local Boltzmann 
weight associated with the Hamiltonian TL in Eq. (1). 11,23 
The DMRG is employed to solve the eigenvalue prob- 
lem 



£T L [aV] * l [<t] = \ l {T)* l [o'] 



(3) 



with X L (T) is the largest eigenvalue of the transfer ma- 
trix and ^ l [ct] the corresponding eigenvector. We em- 
ploy two different boundary conditions: the parallel ones 
(cr 1 = a[ =| and a L = a' L =|) and the antiparallcl ones 
(<r 1 = a[ =t and o L = a' L =|), respectively, for which 
we calculate the largest eigenvalues AjJ(T) and A^(T). 
For the visualization of the spin modulation, we calculate 
the local magnetization 

y\„, *r Ma^r [a] 



Tabic I. Critical temperatures at k = 0.6 known so far. 



Method used 


T 

± c 


r 


Miiller-Hartmann-Zittartz 9 


1.09 




Phenomenological renorm. 8 


1.05 


1.35 


Saqi and McKenzie 19 


1.05 


1.40 


Cluster variation method 20 


0.91 


1.64 


Cluster heat bath method 15 


0.91 


1.16 


Free-fermion approximation 7 


0.907 


1.20 


Non-equilibrium Relaxation 10 


0.89(2) 


0.895(25) 


DMRG (this work) 


0.907 
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Fig. 2. The local magnetization (crj) calculated for L = 118 (i = 
1, 2, ... , 118) with parallel and antiparallcl boundary conditions. 



as a function of position i during the last sweep in the 
zipping process of the finite-system DMRG. 16 We keep 
at most Tfi == 70 block-spin states and vary the lattice 
size from L — 38 to L — 118. Note that under these con- 
ditions the density matrix truncation error 16-18 is kept 
within e < 10~ 13 . 

We use dimensionless units fee = Ji = 1 throughout 
this article. We focus on analysis of the model at k = 0.6, 
where the competing interaction plays an important role 
on the spin modulation. Figure 2 shows the local mag- 
netization at k — 0.6 under and over a transition 
temperature T c ss 0.91 which we will determine more 
precisely. The complete antiphase structure {] ][[} is ob- 
served at T = 0.88 if the parallel boundary conditions are 
imposed (the uppermost) and a twisted pattern created 
by a running domain wall is observed for the antipar- 
allcl conditions (the second from top). The remaining 
two panels display (erj) at T = 0.93, where a modulated 
structure is present for the parallel conditions (the third 
panel) and the antiparallel ones (the fourth). Note that 
the modulation period depends on the applied boundary 
conditions. 

3. Modulation Period 

For the purpose of characterizing the spin modulation, 
we introduce the "domain-wall free energy" 21 



Tbw{T,L) = (-l) n(L h B Tln 



A[ T (T) 



(5) 



where n(L) = [L/2 + 2] + L + 1 represents the 4-site 
periodicity in the antiphase. The Jdw (T, L) represents 
the sensitivity of the free energy per lattice row to the 
boundary conditions. In the antiphase region, !Fuw(T, L) 
exhibits the L dependence 

F DW (T,L)~F DW (T 7 ^) + c(T)L- 2 , (6) 
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Fig. 3. The domain-wall energy /bw(T, °o) at « = 0.6. 



Fig. 5. Behavior of the wave number [fc(T)] 2 = [q(T) - tt/2] z 
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Fig. 4. The domain-wall free energy with respect to size L at T = 
1.45 and k = 0.6 (the symbols). The saw-like fitting functions 
are given by Eq. (8). 



where c(T) is a parameter, which is related to the 'mass' 
of the moving domain wall, and where FnwiT, oo) repre- 
sent the 'stationary' domain-wall energy. Figure 3 shows 
the dependence of !Fr>w(T, oo) with respect to T. The 
domain- wall free energy vanishes at critical temperature 
Tc A ^ — 0.907; the superscript (A) stands for the fact that 
we have estimated the temperature inside the antiphase 
region. 

A detailed analysis of the Tuw(T,L) is required in 
the higher temperature region (T > T c ). The typical In- 
dependence of FuwiT, L) is depicted in Fig. 4. The saw- 
like structure in Fnw(T 7 L) is naturally explained from 
the fact that the system prefers to have 'a natural' wave 
number q(T) for the spin modulation if the system size 
is infinitely large (L — > oo). When L is finite, the bound- 
ary conditions force the system to have a modified wave 
number q'(T), which is quantized as 

q' n (T) = 2mn/(L-£), 

q' n (T) = 2n(m+l)/(L-£), (7) 

respectively, for the parallel and the antiparallcl condi- 
tions, where m is an appropriate integer and £ is an offset. 



As a consequence of the 'forced' shift in the wave num- 
ber, a small increase of the free energy per site occurs and 
is proportional to [q(T) - q' u {T)f and [q(T) - q' n {T)f 
if higher-order corrections are omitted. Paying atten- 
tion to the quantization condition in Eq. (7) and sub- 
tracting T U (T,L) = -fc B TlnA[ T (T) from F n (T,L) = 
— k^Th\ (T), we obtain the saw-like dependence in 
Fdw(T,L) with respect to L shown in Fig. 4. For the 
quantitative determination of q(T), we employ a fitting 
function of the form 

ae~ dL 

^ DW (T,L) = — — {\cos(kL + tp)\ - \on(kL + ip)\} 

(8) 

that contains four temperature dependent parameters: a 
is an amplitude, ip is the phase offset which is related to 
I in Eq. (7), d is a dumping, and 



k(T) = q(T) - ir/2 



(9) 



represents a change of the wave number q(T) from the 
antiphase wave number n/2. 22 (Precisely speaking, as 
shown in Fig. 4, there is a kind of even-odd oscillation, 
and we have to shift <p in Eq. (8) by 7r/4 when L is 
odd.) In Fig. 5, we plot [k(T)] as a function of tem- 
perature T. Performing the extrapolation for k(T), we 
obtain the critical temperature T c = 0.907, which is 
determined in the paramagnetic region. This result is in 
accordance with the previously obtained X"i A ^ • It should 
be noted that the linearity of [&(T)] with respect to T 
is in accordance with the free-fermionic picture on the 
spontaneously created domain walls. 7,8 

Let us observe the temperature dependence of the de- 
cay factor d{T). Figure 6 shows d(T) in logarithmic scale. 
Since our survey is limited to L < 118, the estimated 
d{T) for each temperature contains relatively large fit- 
ting error, which is visible as a fluctuation of the plotted 
data. Among the trial functions we have considered, the 
one 

d(T) = aexp[-p(T-T c )- 1 } (10) 

shows the best fit when a = 1.32(30), f3 = 2.17(32), 
and T c = 0.907(50). Though we don't have any clear 
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Fig. 6. Temperature dependence of the dumping factor. 

picture about the reason of the temperature dependence, 
the obtained T c is in accordance with Tc and Tc . 

As an complementary approach, we calculate 
^bw(T, L) as a function of T at k = 0.6 (see Fig. 7.) 
For simplicity we choose even L in this calculation. 
Again, we obtain the saw- like dependence of Ty>w(T, L) 
which rapidly decays in higher temperature region. It 
is expected that the fitting function in Eq. (8) with 
the previously determined parameters a, k, ip, and d 
explains the plotted data. The curve shown in Fig. 7 is 
thus drawn, and a good agreement is found between the 
plotted data and the fitting by Eq. (8). 

Now, let us observe the 'first zero-crossing tempera- 
ture' of TuwiT, L) from the low-temperature side, when 
L is fixed. Speaking phenomenologically, the wave num- 
ber in the thermodynamic limit q(T) starts to deviate 
from tt/2 at the critical temperature T c . Therefore for 
the finite size system the domain wall energy !Fbw(T, L) 
becomes zero when 

[q(T) - ir/2]L = k(T)L = tt/2 (11) 

is satisfied. Thus the zero-crossing temperature, say T — 
T (L), is slightly larger than T c . From the analysis we 
have performed, we already know that k(T) is propor- 
tional to \JT — T c . Combined with Eq.(ll), we get the 
relation 

^T~rjL)^L-\ (12) 

In order to confirm this, we plot the relative temperature 
r (L)-T (oo) 

T (oo) (13) 

with respect to L~ 2 in Fig. 8, where T (oo) for each k 
is appropriately chosen so that the best linearity in the 
plotted data is realized. As a result, we obtain T (oo) = 
0.907 at k — 0.6, which is in accordance with previously 
obtained transition temperatures. 

For other values of k we obtain T = 1.335 at n = 0.8, 
and T = 1.654 at k = 1.0. The transition temperatures 
thus obtained are used for drawing the phase diagram 
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Fig. 7. Dependence of /bw(^ L) on temperature T at k = 
0.6. The symbols with dotted lines represent the calculated 
•7"bw (T,L) and the full curves are drawn by Eq. (8) . 




L 2 

Fig. 8. The relative temperature t with respect to L~ 2 for k = 
0.6, K = 0.8, and k = 1.0. 

shown in Fig. 1. It should be noted that the divergence 
of the correlation length at the transition temperature 
leads to the relation t oc L~ x l v , 23 and from the relation 
we obtain a critical exponent v = 0.5. This value is in 
accordance with the free-fermion picture. 7, 24 

4. Conclusions and Discussions 

In conclusion, we have applied DMRG to the 2D 
ANNNI model. Observing the spin modulation period, 
we confirm that the free-fermion picture well describes 
the phase transition from the antiphase to the modu- 
lated state. 

Since the correlation length is far longer than the size 
of systems in our study when the temperature is slightly 
higher than T c , we cannot directly judge whether there is 
a 'stable' IC phase or the modulated states are actually 
decaying in the long distance limit. At least we can say 
that we observe no conspicuous singularity in the mod- 
ulation period above T c as shown in Fig. 5. Comparing 
the fact with the possible temperature region for the IC 
phase reported by Shirahata and Nakamura, 10 we con- 
jecture that there is no IC phase in 2D ANNNI model, 
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which is one of the possibility pointed by Shirahata and 
Nakamura. 
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